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1 Introduction 



Entanglement is a fundamental property of quantum systems. Historically, there has been great 
interest in developing efficient theoretical measures of entanglement, one of which is known 
as the bi-partite entanglement entropy [1]. This entropy measures the amount of quantum 
entanglement, in some pure quantum state, between the degrees of freedom associated to two 
sets of independent observables whose union is complete on the Hilbert space. 

To make this definition more precise, let us consider a one-dimensional quantum spin chain 
such as the one depicted in Fig. [T] and let |gs) be the ground state of that chain. We will 
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Figure 1 : A region A of length L of a quantum spin chain 



subdivide the chain in two regions, A and its complement A. In the figure, region A is the 
block of L spins in red whereas A is the rest of the chain, in blue. Associated to this quantum 
spin chain there is a Hilbert space which can be expressed as a tensor product of local Hilbert 
spaces associated to its sites. This can be written as a tensor product of the two Hilbert spaces 
associated to the regions A and A: 

n = A^A. (1.1) 

The bi-partite entanglement entropy is the von Neumann entropy of the reduced density matrix 
Pj^ associated to A, that is: 

SA = -TrA{PA^ogpA) , p^ = Tr^(|gs)(gs|). (1.2) 

It is the entropy of the region A with respect to the rest of the chain, regarded as its environment. 
Another measure of entanglement that is widely studied in the literature is the Renyi entropy 
[2], given by 

whose n — 7- 1 limit gives the von Neumann entropy. 

The bi-partite entanglement entropy of one-dimensional (integrable) quantum spin chains has 
been extensively studied in the literature from different points of view. One may for instance 
consider a quantum spin chain divided in two parts, both of which are infinitely long. If the spin 
chain has a finite gap, the entropy saturates to a constant value which depends on the parameters 
of the model under consideration and can in some cases be computed analytically [31 [U [5]. 
The methods used are based on the algebraic Bethe ansatz or on the corner transfer matrix. 
Another interesting set-up are quantum spin chains which are infinitely long but where the block 
A is kept of finite length. In the critical regime, it has been found through a combination of 
numerical and analytical approaches (many of which can only be applied for integrable models) 
[3 El El ini [HI [E] that the entropy of a block diverges logarithmically with the size of the block, 
as expected from conformal field theory (CFT) considerations \13 \ I14 ^ \TE \ I16j. 

The regions A and A can take many different shapes and need not necessarily consist of 
connected blocks of spins. That is, we may also consider a situation in which both A and 
A are made out of the union of several disjoint blocks of spins. The entropy of disconnected 
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regions has been the object of numerical and analytical study recently in the context of CFT 
[ni dSl da EOl EI] and quantum spin chains and lattice models [22l [23] . 

The goal of this paper is two-fold. First, we introduce the quantum chain precursors of 
branch-point twist fields of quantum field theory (QFT), and explain their use for the computa- 
tion of the bi-partite entanglement entropy. Branch-point twist fields were introduced in a series 
of recent works in order to compute the bi-partite entanglement entropy of 1+1-dimensional 
QFT [MlEHlEelETlEH]. Using the "replica trick," one writes the entropy ([Oj) as 

SA = -lhnJ-TrA{p'X). (1.4) 

The quantity Tr^ (Pa)) for integer n, can be interpreted as the partition function associated to n 
copies or replicas of the model, branched in a particular way. In the context of l+l-dimensional 
QFT, we found [23] that the partition function is proportional to the two-point function of a 
pair of branch-point twist fields T, T: 

Tr^(PA)~(r(r)r(0)). (1.5) 

These twist fields are characterized by their exchange relations with any local field </3j of the i-th 
copy of the mode0 

^MTix) = Qix' > y') T{x)^,+M + e{x^ < y') T{x)^Ay) (1-6) 
ipMf{x) = e{x' > y') f{x)^,^M + 0(^' < y^) f{x)^i{y) (1.7) 

for i = 1, . . . , n and where we identify the indices n + i = i. Hence, branch-point twist fields are 
local twist fields in the replica model associated with generating elements of its Z„ symmetry. 
Locality plays a crucial role in evaluating their two-point function. Of course, in order to obtain 
the entanglement entropy, one performs a continuation to real n and the limit n — )■ 1 in the 
resulting expression. 

In the context of quantum spin chains, it is possible to define cyclic permutation operators 
in terms of which the equivalent of the branch-point twist fields of QFT can be constructed. 
Correlation functions of such permutation operators yield the entanglement entropy of subsets 
of the chain. In this context, there is also a concept of locality associated to such permutation 
operators, although we will not make explicit use of it. We will rather show, through examples, 
that the use of cyclic permutation operators provides a clear combinatoric method for evaluating 
the entanglement entropy. 

Second, we use this tool in order to compute the von Neumann and Renyi entanglement 
entropies of the state obtained in the limit A — t- — l"*" of the infinite-length anisotropic XXZ 
quantum spin chain (A being the anisotropy parameter). The limit A — t- — 1+ and some special 
entanglement features of the associated ground state have been previously considered in [6], 
where the concurrence was computed, albeit for the case of finite chains. In the case of infinite 
chains, there are subtleties involved in this limit. In order to properly address them, we provide 
a treatment of quantum states in infinite-length chains based on linear functionals on the space 
of finitely-supported operators. Using a completely natural compact convergence topology on 
linear functionals, we describe the limit A —t- —1"*" as a proper quantum state in this framework. 

We interpret our results by writing the state as a "hybrid" between a factorisable ground 
state of the XXZ chain at A = —1, and a state preserving the antiferromagnetic property of 

*Here we employ the standard notation in Minkowski space-time: x'^ with u = 0,1, with being the time 
coordinate and being the position coordinate. 
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XXZ ground states in the range —1 < A < 1 (the z component of the total spin being zero). 
We propose the existence of a new length scale, below which the spins behave as if they were in 
a ground state of the XXZ chain at A = —1, and above which this antiferromagnetic property 
holds. As A — )• —1''', this scale tends to infinity, and the entanglement entropy of finite blocks 
only measures the entanglement due to the antiferromagnetic condition. 

Recent results [5] show that the entanglement entropy between the semi-infinite halves of 
the infinite-length XYZ spin chain behaves, near to certain critical lines, in ways that are not 
explained by the usual CFT arguments. Although the cases that we study do not overlap with the 
cases studied there, we obtain similar conclusions. The von Neumann and Renyi entanglement 
entropies associated to finite blocks of sites in the A — )• — l"*" limit-state behave, as the blocks 
become large, in a way that agrees with scale invariance, but not with conformal invariance. We 
suggest that some results of [5] may be due to an interplay between the new length scale that we 
propose and the usual length scale associated to the mass gap. We hope that the relationship 
between our results and those of [5] will be made more precise in a future publication |29j . 

The paper is organized as follows: in Section 2 we introduce a new type of local operators 
which can be defined for a replica version of any quantum spin chain model. These are local 
cyclic replica permutation operators, and we provide an explicit expression for them in terms 
of elementary matrices. For completeness, we provide proofs of their main properties, including 
their action on quantum states and their exchange relations with respect to other local operators 
in the chain. We explain how the bi-partite entanglement entropy can be expressed in terms of 
correlation functions of such operators. In Section 3 we apply the results of Section 2 to the 
particular case of the spin-^ XXZ quantum spin chain and compute the von Neumann entropy 
associated to one or two spins in the infinite chain. In Section 4 we define the notion of quantum 
state in infinite-length spin chains. In Section 5 we find exact formulae for the von Neumann and 
Renyi entropies associated to the quantum state occurring in the limit A — t- — l"*". We study the 
asymptotic behavior of the entropy for large subsystem size and find that it scales logarithmically 
but not conformally. We finish this section with a discussion and some speculation about the 
nature of this quantum state. We conclude the paper in Section 6. 

2 Replica permutation operators and entanglement entropy in 
quantum chains 

A quantum chain is a quantum system physically composed of a number N of locally-interacting 
sites along a one-dimensional chain. Mathematically, two conditions are imposed. First, the 
Hilbert space is Hat := ®f=]Vi where Vi = represents the states of the individual site i, 
with the usual inner product (we may think of this as a spin chain of spin {d — l)/2j3. Second, 
the Hamiltonian H = Yl^i ^i' representing the interaction between the sites, is a sum over 
local energy operators hi, usually without explicit dependence. Each hi factorises to the 
identity on all but a certain number (independent of A^) of sites. This gives rise to a notion 
of "neighbourhood": a site j is in a neighbourhood of a site i if hi does not factorise to the 
identity on site j. It is this notion of neighbourhood that tells us that the chain is actually 
one-dimensional, that provides a concept of distance (e.g. the minimal length of a sequence of 
sequentially intersecting neighbourhoods required to cover two given sites), that tells us about 
its most natural topology (e.g. a chain with boundaries, or a periodic chain), and that usually 
gives rise, in the large-A^ limit, to factorisation of correlation functions at large distances. This 

^It is not important for many of the considerations that all sites have the same dimension; yet, for simphcity, 
in most of this paper we will specialise to ci = 2. 
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notion of neighbourhood is also at the source of the emergence of a local QFT in the scaling 
limit, when a quantum critical point exists. The notion of neighourhood extracted from the 
energy density can often be made to agree with the usual notion of neighbourhoods on a chain 
(and we will assume this throughout). 

In this perspective, it is natural to look for a quantum-chain precursors for our branch-point 
twist field defined in 1-1-1 dimensional QFT [23]. Recalling the exchange relations ()1.6p - ()1.7p . it 
is clear that these will be certain permutation operators (elements of a linear representation of 
the permutation group), and like in the context of QFT, the entanglement entropy in quantum 
spin chains can be expressed in terms of correlation functions of these operators. 

We consider n independent copies, or replicas, of a quantum chain. The Hilbert space can 
be described by 71%'' ^ TiP := 0^=1 Va,i = where Va,i = and the vector 

(n) (n) 

space at site i is := 'i^a=i^oi,i- It will be convenient to denote vectors belonging to 

by |siS2 • • • Sn)i, Sa G {1,2,..., d}, where the "spin" Sa, a = 1, . . . ,n belong to copy a of the 

model at that particular site i. 

The precise permutation operators that we need, denoted %, are local cyclic replica permu- 
tation operators, acting non-trivially only on site i in the n-copy spin chain, with i = 1, . . . ,N. 
They act as follows: 

Ti\siS2 ■ ■ ■ Sn)i = |S2S3 . . . S„Si)j. (2.8) 

For A C {1, 2, . . . , A^} a set of sites on the chain, we will also use the notation 

Ta:=\{%. (2.9) 

Following similar ideas as in QFT, the trace of the n-th power of the reduced density matrix 
associated to a state IV') ^T-Ln is 

{|<j>f >e8} 
1 " 



(2.10) 



where |^') = IV')'^") and the sums are over orthonormal bases. Hence, the bi-partite Renyi and 
von Neumann entanglement entropies of the block A can be evaluated as follows: 

and 



= - lin, A C-™) . (2.12) 

(Naturally, any other cyclic element of the permutation group could have been used, by permu- 
tation invariance). As usual, an analytic continuation in n is understood in the latter expression. 
We will come back below to the full relationship between these permutation operators and QFT 
branch-point twist fields. 
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2.1 Expressions in terms of elementary matrices 

We now study the relation between these permutation operators and elementary matrices. This 
will be useful because correlation functions of the latter have tractable expressions in the context 
of integrable quantum chains. For simplicity, let us consider the case of a spin-| chains, where 
the spin variables si, . . . , Sn above can only take two values (that is, d — 2). The generalisation 
to higher spins is straightforward and explained below. The 2x2 elementary matrices , for 
e, e' S {1,2}, have matrix elements given by 

(-£'" )kk' = Se,kSe',k'- (2.13) 
From these matrices, we form the operators 

E<, a = l,...,n and i = l,...,N, (2.14) 

which act on site i and copy a of the quantum chain as the matrices E''^' , and everywhere else 
as the identity operator. 

In the simplest non-trivial situation, the two-copy model (n = 2), the permutation operator 
Ti simply exchanges copies 1 and 2 at site i. It has a well-known expression in terms of elementary 
matrices given by 

T = E\\El\ + E\lEl\ + El}Ell + El}Ell (2.15) 

It is a rather easy (although tedious) exercise to generalize "by hand" this result to higher values 
of n, keeping spin-^. We show here just the results for n = 3 and n = 4, 

Z^ll Z^ll I ZT'l^ T7'21 1722 I T7i21 T71I2 TTill I I7'22 T7i22 T7i22 i 

i - ^l,i^2,i^3,i + ^l,i^2,i^3,i + ^l,i^2,i^3,i + ^l,i^2,i^3,i + 

TTill 7721 Tpl2 I Tpl2 TTill Tp21 , T7'22 Tpl2 Tp21 1 Tp21 Tp22 Tpl2 /r) -1 \ 

^l,i-^2,i^3,i + -^^1 ^2,i^3,i + ^l,i^2,i^3,i + ^l,i^2,i^3,i^ l^-J-Oj 



-r — f?ll f?ll _L f?21 pl2 Tpll Tpll I TpU p21 77112 pll I 77121 77n22 77nl2 pll , 

TT'll iT'll I7'21 Tpl2 I Tp21 Tp22 Tp22 Tpl2 I 77'21 rpl2 Tr21 Tpl2 , rpll Tp2\ Tp22 Tpl2 i 
^l,i^2,i^3,i^i,i + ^l,i^2,i^3,i^i,i + ^l,i^2,i^3,i^i,i + ^l,i^2,i^3,i^i,i + 

Tpl2 jpll TTill Tp21 I T7'22 Tpl2 jpll Tp21 1 77'12 77'21 Trl2 Tp21 I rr'22 rr'22 Tpl2 Tp21 . 
^l,i^2,i^3,i^4,i + ^l,iJ^2,i^3,i^4,i + ^l,i^2,i^3,i^4,i + -C'l,i-C'2,j-C/3,j-a4,j + 

T7'12 Tp21 Tp22 Tp22 I 77'22 Tpl2 Tp21 Tp22 I 77'12 TT'll 77'21 r7i22 I 77'22 7T'22 77'22 7T'22 /r) -1 '~r\ 

^l,i^2,i^3,i^4,i + ^l,i^2,i^3,i^4:,i + ^l,i^2,i^3,i^4,i + ^l,i^2,i^3,i^i,i- l^- ' j 

Once a few explicit cases have been worked out the general structure of the permutation operator 
for generic n quickly starts to emerge. There are several features that we could have predicted 
from the start. For example, all the formulae above are sums of exactly 2" terms, which is also 
the number of distinct basis vectors |siS2 ■ ■ ■ Sn)i that can be constructed for a spin-| model. 
The permutation operator is a sum of 2" such terms since its action ()2.8p is a one-to-one map 
between such basis vectors. In addition, all expressions are symmetric under the combined 
exchanges E^^ o £'^^ and E^"^ ^ E'^^ . 

One can characterize the precise 2" terms that will appear in the expression for the permu- 
tation operator through a set of four simple rules: 

Rules 2.1 The operator % is the sum of each possible term, with coefficient 1, that is a product 
of matrices E^f"', a = 1, . . . , n respecting the following rules: 

1. a matrix E^- can only be followed in the product by E^_^-^ ■ or E'^^_^^ ■, 

2. a matrix E'^- can only be followed in the product by -E^+i j or E^_^^ ■, 
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3. a matrix E^- can only be followed in the product by E'^_^-^ ■ or E^^^-^ ■, 

4. a matrix E'^\ can only be followed in the product by -E^+i j or E'^_^_^ ■ 

(with n + 1 = 1 and where the rules apply cyclically in the product: the first factor follows the 
last factor). 

This restricts greatly the type of terms that will emerge. For example there will never be terms 
which involve both matrices E^^ and i?^^ and no other type, although the terms £'^^£'2^ • • • E^\ 
and £'^^£'1^ • • • E^^- are allowed and always appear, as we can see for n = 2, 3, 4 above. Similarly, 
the number of matrices £^^ and £^^ in a given product is always the same. In fact, the rules 
prescribe the following simple structure for the elementary matrices at any given site along the 
n copies: strings of any number (including zero) of E^^ and £^^ are separated by £^^ (in the 
junctions from £^^ to £^^) and by £^^ (in the junctions from £^^ to £^^), with a condition of 
periodicity. This point of view will be very useful in explicit calculations below. 

The rules above immediately give an expression for the permutation operator of the form 

2 

ei,...£n = l 

This reveals the structure of the trace of some matrix. More precisely, if we introduce an 
auxiliary space 14ux — as well as the operators 

7«,i;aux = ^ ^lux^a^i = ( j^l2 j^22 ] ) (2-19) 



=1 



which acts non-trivially both on the auxiliary space, and at site i, copy a, then the permutation 
operator at site i can be written as the following trace 

/ n \ 

"Ti — Tr^ux (^l,i;aux^2,i;aux ' ' ' ^n,i;aux) — Tl^ux I ^o,i;aux 1 (2.20) 

\a=l J 

where the product is, by convention, ordered from left to right in order of increasing a. It is 
easy to obtain (|2.15p . (j2.16p and (j2.17p from (|2.20p by setting n = 2,3 and 4, respectively. As 
it should be, Ti reduces to the identity matrix for n = 1. 

We note that the operator Ta^i-aux is simply a permutation operator exchanging the auxiliary 
space with the space (a, i). The form (j2.20p of the cyclic permutation operator is something that 
is well-known in the context of integrable quantum spin chains solved by means of algebraic Bethe 
ansatz. In this context, the £-matrix is, at a certain value of the spectral parameter, simply 
a permutation operator exchanging an auxiliary space and the space associated to a given site 
of the quantum chain (see, e.g., [30]). The monodromy matrix is the trace of the product of 
L-matrices, exactly of the form (|2.20p . and is identified, at this special value of the spectral 
parameter, with the translation operator along the periodic chain. In our context, however, the 
operator that we obtain permutes the copies of a replica model at a given site, rather than the 
sites of the chain. 

In fact, this connection between our expression ()2.20p and the algebraic Bethe ansatz becomes 
very explicit if one employs the solutions to the inverse scattering problem found in |31j to express 
the operators E^- in terms of the entries of the monodromy matrix. We suspect that this 
relationship could be useful in the future for studying the properties of the replica permutation 
operator within the Bethe ansatz framework. 
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2.2 Exchange relations 



In the derivation above, it is clear that the Rules 12.11 the expression (j2.18p . and the trace 
expression (|2.2U|) are equivalent. A careful analysis of the permutation action (j2.8|) would also 
give rise to the four rules stated. Here, for completeness, we will show explicitly that equation 
(j2.8|) is a consequence of (|2.2U|) . In order to do so, we will show that the main exchange 
property of replica permutation operators (the property that constitutes the starting point for 
the definition of twist fields in QFT) is satisfied by the operator (j2.20p . 

Given a 2 by 2 matrix O, let Oa^i be the operator on the quantum chain that acts as O on 
site i, copy a, and as the identity operator everywhere else. Then, we have the following: 



Lemma 1 The operators % defined by \2.2U\) satisfy the relations 
(with Oo,i = On,i) for all i = 1, . . . , N , all a = 1, . . . ,n, and all O. 



(2.21) 



Proof. This relation follows from the trace expression ()2.20|) along with two identities. Denoting 
by Oaux the operator acting non-trivially as O on the auxiliary space and as the identity operator 
everywhere else, we only have to show that 



On 



O T 



a,2;aux>-^aux 



Oa 



O T 



The first identity can be derived as follows. With 



(2.22) 



O 



we have 



T O 



?2l 



Ea,i E'^,i 



nl2 



aux 



On oi2 

021 022 



Oa,^ 
Oa,i 
■121 



02lE^\ + 022-Ei^i 02lE'f^j + 022-E, 



aux 

oiiE-i], + ouEf. 

21 I ^ i?22 



Ei]iOa,i 



which is equal to 



o.. 



T ■ 

aux Q,j;aux 



The second relation of ()2.22p follows in a similar way. 

The defining property (j2.8p is essentially a consequence of (j2.2ip . 



El]iOa,i 
22 Q . 



e: 



(on o^2\ ( E^, E^, \ 

02 J. Am <Jaux 

/ ouEi], + o,2El], oiiEl], + ouE^l \ 

V 02lEi], + 022Ei'^, 02lEl]^ + 022Ef.^ ) 



(2.23) 



(2.24) 



(2.25) 



Lemma 2 The operators % defined \2. 20]) satisfy the relations 

Proof. We start by writing the quantum state of the spin-^, n-copy quantum spin chain at site 
i in terms of elementary matrices -E^^^^. Let |0)j be a reference state at site i for which all spins 
of all n-copies are up (si = • • • = =t)- Any other of the 2" possible spin configurations can 
be generated by acting on this reference state with a combination of elementary matrices. For 
example, if we want to lower the spin of copy 1, we just have to compute £^^^JO)i, whereas acting 
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with EI\ will leave the state unchanged and the action of Ef^ and E'J^ gives zero. Therefore, a 
generic state is written as 



\SlS2 



■Sn)i= \YlE'jf]\0)i„ js^ = l for sa=t and j,^ = 2 for s„ =i, (2.26) 



Ka=l 



and from the expression (|2.18|) . we have 7i\0) = |0) since the only non-zero summand is the one 
where = 1 for all a. We now evaluate the action of the permutation operator on this state 
by successively employing ()2.2ip . 



n 



ri\si,S2,...,Sn)i = % (^n^fcj 

= ( n ^a-M ) "7^10)' = ( n ^"-M 1 = 1^2, ■ ■ ■ , ^n, 5i>i. (2.27) 



This establishes (12. 8D. 



>^a=l / \ci=l 



2.3 The conjugate permutation operator T 

Other permutation operators, % at each site i, take us from copy a to copy a — 1. We can 
simply obtain 71 by reversing the order of the replicas, that is 

Ti — Tr^ux (^r!,,i;aux^n— l,i;aux " " " ^l,i;aux) j (2.28) 

or, alternatively, taking the trace, 

2 

fi= V E^^'^^E^^^'i ■E^^% ■■■■ (2.29) 

ei,...e„=l 

One may check this explicitly by showing that 71% = 1. This can be done by using the following 
property of the elementary matrices: 

^^f^sr = ^..,e3£^^r- (2.30) 

Employing the representations (j2.18|) and (|2.29|) for the permutation operators (including some 
extra matrices in the products for clarity) we have 

2 

lili - 2^ ^ -^„-l,i^n-2,i • • • -^3,i ^2,1 El,i 

ei,...e„ = l 

2 



^ \ ^ J^^2^1 J^^4^3 . . . J^^n.-l^n-2 j^^l^n 

1,2 2,i 3,i Ti — 2,i Ti — 1,2 

f' e' —1 

ei,...e„ = le[,...e'„ = l 
2 

ei,...e„ = l 

The last equality follows from the fact that all matrices E^^ are diagonal, so that each product 
above is a diagonal matrix with a single non- vanishing entry 1. The sum over all such matrices 
(there are exactly 2" of them) is the identity. 
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2.4 Replica permutation operators for higher spin chains 

We can now generalize the results above to more general spin chains, with arbitrary spin repre- 
sentations living at each site of the chain. Indeed the exchange relations (|2.21|) which we have 
just proven can be employed as a starting point for finding a realization of the permutation 
operator which is valid for higher spin chains. It can be shown that, assuming a dimension-^ 
representation lives at site i of the quantum spin chain, the permutation operator 71 can be 
expressed exactly as in (|2.2U|) with Ta^j^aux given by 



'a,j;aux — 





tp22 










■ K', ) 


aux 



(2.32) 



It is not too difficult to show that all properties established in the previous section for spin-^ 
also hold for the general case. We will not dwell on such proofs here as in the remainder of this 
paper we will concentrate on the spin-^ case. 

2.5 Local-unitary operators and twist fields 

Since the operator Ti is real, we immediately find, from (|2.20p and (|2.28p . that 

T — 
ii — li ■ 

Hence, 7i is a unitary operator. It is of course the local element of the unitary operator associated 
with the cyclic replica permutation symmetry of the n-copy Hamiltonian 

N N n 

^("^ = E^!"^ = EE^- (2.33) 

(where h^^i, for a = 1, . . . , n and i = 1, . . . ,N, acts as hi on ^jLi^aj, and as the identity on 
^jLiVa'j \/ a' ^ a). More precisely, the unitary operator that is the product of % over all sites 

i, i.e. T{i^,,,^N} = rii^i 'Ti, implements the cyclic replica permutation of the chain, and commutes 
with the n-copy Hamiltonian. A unitary operator that can be written as a product of single-site 
unitary operators is sometimes called a local-unitary operator. 

In fact, the local-unitary operator 7{i^,,,^/vf} commutes also with the energy density, 

T[i_N},h^^] =0. (2.34) 

In general, it is expected that any unitary operator associated to an internal symmetry, i.e. that 
is local-unitary, commutes not only with the Hamiltonian, but also with the energy density. To 
such an internal symmetry, one may then associate twist fields; e.g. in the case of 7i, 

N 



'7{i,i+l,...,Ar} — JJ^i- 



A standard example is the twist field involved in the construction of fermion operators on the 
Ising spin chain. Quantum chain twist fields give rise to the usual QFT twist fields in the scaling 
limit. Twist fields have the property of being local "in the bulk": their commutators with the 
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energy density vanish for large enough separations (and away from the end of the chain N), by 
virtue of their association with an internal symmetry transformation. In the case of the replica 
permutation symmetry: 



[Tj, hj] = for j far enough from i and from N. (2.35) 

This has important implications in QFT, and is likely to have important implications as well 
in the context of integrable quantum spin chains. We will not investigate further here such 
implications. 

Consider the permutation operators 71 and their twist fields Tj. In the scaling limit, this of 
course gives rise to branch-point twist fields, associated with branch-point singularities on the 
surface where the QFT is defined. Like in QFT, the entanglement entropy (and likewise the 
Renyi entropy) of a region A of the quantum chain can be expressed via a correlation function 
of as many pairs of twist fields Tj, t| as there are connected sub-regions. For instance, for a 
single connected region A starting at site i and ending at site j, we have 

^Renyi 1 , ( {m^\ 



2.6 Local-unitary invariance 

One of the properties of the entanglement and Renyi entropies is that they are invariant under 
local-unitary transformations (transformation of the state by a local-unitary operator). The 
precursor of this invariance is of course that the replica permutation operators Ta are invariant 
under tensor-powers of local-unitary transformation^, U®"'Ta{U^)^"' = Ta for U a local-unitary 
operator on H. In general, we expect that if two quantum states have the same Renyi entropies, 
then they are related by a local-unitary transformation. 



3 The XXZ Heisenberg spin-^ chain: main features 

In order to illustrate the technique proposed, we now choose the XXZ Heisenberg spin-^ chain 
as a benchmark. It is one of the most studied quantum chains and includes other interesting 
theories as special cases. The XXZ Heisenberg spin-^ finite chain of length N is characterized 
by the Hamiltonian 

N 

Ha = Y. + ^M+i + ^(^^.Vi - 1)) . (3-36) 

where aj'^'^ are the Pauli matrices associated to site j of the chain and acting on the two- 
dimensional space Vj = C^, a spin-1/2 module. We consider periodic boundary conditions 
^AT+i = o'l'^'^, and for simplicity we restrict to being even. The parameter A is known as 
the anisotropy parameter and, depending on its value, the physical properties of the model can 
change dramatically. In this paper, we will be looking at the region — 1 < A < 1. In this region, 
the model is antiferromagnetic and preserves the z-component of the total spin: its ground state 
has total 2;-spin equal to zero. As — t- oo, the model is massless: the gap between the unique 
ground state and the first excited state tends to 0. For A = 1 the model becomes what is 

"^The set of (products of) replica permutation operators associated with all permutation elements spans the 
space of local- unitarily invariant operators (but these operators are not all linearly independent). 
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known as the XXX model or Heisenberg spin chain, whereas for A = it reduces to the XX 
model, which possesses a free fermion description. For A = — 1, the ground state of the model 
is very different. It is unitarily related to the ground state of the ferromagnetic Heisenberg spin 
chain, and is not unique. In all cases, the ground state is invariant under translations (cyclic 
permutations of the sites along the chain). The anisotropy parameter is commonly expressed in 
terms of another parameter rj as 

A = cosh(77), (3.37) 

with < < vr. 

The characterization of the physical states, including the ground state, of this and many other 
integrable quantum spin models can be carried out very effectively by means of the algebraic 
Bethe ansatz approach [32l [30]. The details of the approach can be found in many places 
(a particularly good review can be found in [30]) and we will not review those here. The 
entanglement entropy is a function of the correlation functions of the model. The evaluation 
of correlation functions has been performed in the literature by various methods: using the 
algebraic Bethe ansatz [33l [30l [3l] with the solution to the inverse scattering problem [311 |35l [36] , 
and through a computation of form factors [37\ [38l [39] in quantum spin chains of infinite length 
via g-deformed vertex operators. 

We will denote by |V')a the ground state of the Hamiltonian ()3.36p . We will also use the 
notation E^^ for the operator on T-L acting as the elementary matrix (|2.13|) at site i and as 
the identity elsewhere. Instead of considering correlation functions of Pauli matrices, we will 
consider fundamental blocks, correlation functions of (i.e. quantum averages of products of) 
elementary matrices in the ground state |V')a- We will use the following notation (with implicit 
dependence on A): 

a(vi n Ei-<\i^)A 

= , (3.38) 



(11 



where A C {1, 2, . . . , N} is an index set. 

We will denote by |^') (with implicit dependence on A) the vector in the replica model "H*-""^ 
corresponding to the tensor product |V')a" ground state. This is of course the ground 

state of the n-copy Hamiltonian H^^\ By factorisation, we have, for instance, 

= Ui^r')- (3-39) 

^ ' ' a=l 

We note that it is a simple matter to verify the general locality property ()2.35p of the twist 
fields Tj in the present case. With hj = a^aj^^^ + a-^(T^j_^_i + A{ajaj_^_-^ — 1) and /i^"^ defined as 
in (j2.33p . we find 



J — "i— 1 J ■■■ ) 

a=l 
for j / N. 

3.1 Entanglement entropy of one and two sites 

We now perform simple computations of the entanglement entropy between one or two sites 
and the rest of the chain using the permutation operator technique. This illustrates how the 
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technique directly gives the combinatorics for evaluating the trace of the reduced density matrix 
in terms of correlation functions (avoiding the explicit diagonalisation of the matrix). 

The one-site computation can be thought of as a consistency check of our formalism, since 
the outcome is known a priori. This is because the ground state of the spin-^ XXZ chain is 
characterized by equal numbers of spins up and down. Hence if we pick one spin in the chain, 
its projection in the z-direction is entirely determined by the values of that projection for the 
remaining spins in the chain. In other words, the entanglement between one spin and the rest 
is maximal and it is well known that this maximum value is exactly log(2). 

The entropy that we want to compute is 

S,„ = - ,i,„ ^ (!^) . (3.40) 

{En = {Ef') = 0, (3.41) 



In order to evaluate this, we note that 

^12 \ /i7i21 

and 



{E}') = {Ef) = 1. (3.42) 

The first equality (|3.4ip is due to the fact that the operators Ef^ and Ef^ have the effect of 
reversing the spin at site i, turning the quantum average above into the scalar product of the 
ground state with another orthogonal state. The second equality (j3.42p is due to the fact that 
the quantum averages involved represent the probability of finding the spin at site i up, and 
that of finding it down. These probabilities are exactly ^ due to the nature of the ground state. 

Factorisation ()3.39p and the expression ()2.18p mean that we have to evaluate a sum of 
products of one-point functions. Because of (j3.4ip . the only non- vanishing contributions to the 
quantum average of 71 will come from the two terms that involve only matrices E^^ or i?^^. 
That is: 

i^^^ = {EjY + {EfY = 2'-\ (3.43) 
The entanglement entropy is therefore, as expected 

Su\ = - lim -^2i-" = lim 2^-'^ log(2) = log(2). (3.44) 

We now consider the computation of the entanglement entropy of spins sitting at sites 1 and 
m + 1 with respect to the rest of the chain: 

(this is the same as S'jj „,,4-i} thanks to translation invariance). 

Given the factorization property (|3.39p . we know that 5*11 „i_|_i}. will be given in terms of 

two-point functions involving pairs of operators E^^''^ E^'_^\. However, as before, many of these 
correlation functions are vanishing. More precisely, they do so whenever an operator E^'^ or E'^^ 
appears in combination with E^^ or E'^'^. There are only six non- vanishing two-point functions: 

(-^i^^m+i) = {Ej'^ E^j^^) , {EI^E^_^^) = {Ef EI:^_^_^) , (3.46) 

and 

= (3-47) 
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where the equahties are again due to spin reversal symmetry. This means that although the 
operator TiTm+i is a sum of 4" terms, each of which involving 2n matrices E^^ , inside the 
correlation function most of these terms will give a vanishing contribution. 

We may evaluate the correlation function by using Rules 12. H recalling the structure men- 
tioned there. Here, we may consider a square lattice of cylindrical topology, with 2 columns, n 
rows, and periodicity along the columns. Each cell is occupied by an elementary matrix: this 
represents a term in the product 7iTm+i- Rules [27T] imply that in order to determine a column, 
we only need to provide the positions of matrices E^^ and E'^^, with the condition of alternation, 
unless there are none. In the latter case, there are only two possibilities: whether the column is 
filled with E^^, or it is filled with i?^^. On the other hand, the constraints about non- vanishing 
correlation functions imply that on a row, if a E^"^ is present then a E'^^ must also be present. 
This implies that in order to determine a term, we only need to determine the positions of 
matrices E^"^ and E"^^ along one column only - say the first column - unless there are none. 

Suppose that along the first column there are s matrices E'^'^, and q pairs of matrices E^'^ 
and E'^^. Then there are n — s — 2q matrices E^^. In the case where q = 0, hy the discussion 
above not all values of s are available, and there are exactly four terms: 

For g > 0, all values q from 1 to [n/2] (where [•] means integer part) occur, and all values of s 
from to n — 2(7 occur. By factorisation, every term characterised by such s and q gives rise to 
the same product of correlation functions: 

/pll 17.22 \n-s-2(j/ p22 pll \s/pl2p21 \g/p21pl2 \q 

In order to count the number of such terms, consider two possibilities: whether the jump from 
n to 1 along the column breaks a string of (possibly a vanishing number of) matrices i?^^, or it 
breaks a string of matrices E^^. In the first case, the number of terms is pq+i{s)pq{n — s — 2q), 
and in the second, it is pg+i(n — s — 2q)pq{s), where Pq{s) represents the number of partitions 
of s into a sum of q parts, with the number is also included as a possible part. Hence, the 
number of such terms is 

Cn{s,q) = pq+i{s)pq{n - s - 2q) +pq+i(n - s- 2q)pq{s). (3.48) 

A formula for Pq{s) can be easily obtained by noticing that the function (1 — x)"" is precisely 
the generating function of such coefficients, namely 



^ oo 



/ 1 \ 



where ^ ^ ^ ~ is the binomial coefficient. We can easily show that the coefficient ()3.48p 

can be written as 



n — s — q \ 5—1 

Employing this result and from the discussion above, we find that the correlation function that 
we want to compute is 

^IT^^^ ^ = Y.Y.Cn{,,s){El^E^^^,r-^-^'^{EfE^^,Y{E\^E^^^,Y 

^ I ' q=l s=0 

+{Ei'E'j^,r + {E?E'j^,r + {Ei'E'j^.r + «i?^Vi>"- (3.51) 
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The sum above is relatively complicated to compute. However, it simplifies greatly if we 
employ the equalities ()3.46p - (|3.47p . We then obtain 



{^\TlTm+l\^) 



n—2q 



q=l s=0 



It is now possible to carry out the sum. The sum in s yields simply 



n—2q 
s=0 

Finally, the sum in q gives 

/Eillp22 \n-2g/pl2p21 \2q 



n 
2q 



(3.52) 



(3.53) 



q=l 



n 
2q 



= -2{E\'E^_,,r + {{El'Ei\,) + «i?^Vi))" + {{El'Ef:^^,) - (i^^i^^V,))" (3.54) 
Therefore, the final expression for the two-point function of permutation operators is 

^^'^Jl^y'"^^ = m'E^^^,) + (i^Pi^^Vi))" + m'E^^^,) - {El^E^J^^,)r + 2{El^E^^^,r, 

(3.55) 

and the entropy becomes 



-((-^l^^m+l) - (-E'P-E'm+l)) log [(-E'l^-E'm+l) ~ 

-2{El'E'Jl^,)log{El'E},l^,). 



(3.56) 



Note that in (|3.55p we can immediately read-off the eigenvalues of the reduced density ma- 
trix (the three quantities that are taken to the power n) and their degeneracies (the integer 
coefficients of these powers). 

It is simple to rewrite this expression in terms of Pauli matrices. Employing the short-hand 
notation (cfCTm+i) — z{m) and {(ytcy^j^i) = s{m), we have 



(1 — z{m) + 4s(m)) 



(1 



z[m 



4s (m)) 



log 



log 



[1 — z{m) + 4s(m)) 



(1 — z{m) — 4s(m)) 



(1 + z{m)) 



log 



1 + z{m) 



(3.57) 



Prom here let us specialise to the case of infinite length, N ^ oo. This makes sense because 
correlation functions of local operators have finite limits. Then, one feature of the expression 
above is that for m — )• oo the correlation functions s{m), z{m) — )• by large-distance factorisa- 
tion, and therefore the entropy at large distances saturates to its maximum value 



lim 5'{i,.m+i} = 21og(2), 

<y-i — i^r\r~: ^ ' J 



(3.58) 
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which is exactly twice the entanglement entropy of a single spin computed above, as expected. 
One can now evaluate the entanglement entropy of two spins for a great variety of spin-^ models 
for which the correlation functions are known. 

We will conclude this subsection by providing some numerical results for the entanglement 
entropy S^i ^^i^ at particular values of m and A. Although the general expressions for the 
two-point functions involved in (|3.57p are given in general in terms of complicated multiple 
integral representations [391 EHl EH HO] , those integrals have been done at least for small values 
of m and different values of the anisotropy parameter. Very useful tables listing explicit values 
of the two-point functions for many values of A and m = 1,2 and 3 can be found in [4H I42j . 
Employing those results, it is possible to evaluate the entropy for various values of the anisotropy 
parameter as shown in Fig. [2]^a). In addition, the two-point functions of spin operators admit a 
simple form for special values of the anisotropy parameter. One such value is A = or = ^ 
which is commonly known as the free Fermion point (see Fig. El^b)). 
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Figure 2: (a) The entanglement entropy of two spins S^i ,^j^iy in the XXZ chain for different 
values of the anisotropy parameter and m = 1, 2, 3. We observe how the entropy tends to grow 
with m until reaching the value 21og(2) for very large m. An exception to this is the case 
T] = IT. In this case the entropy takes the value | log(2) irrespectively of m. This is due to 
the very special nature of this point, for which the two point function s(m) oscillates between 
the values as m — )• oo (rather than vanishing). The limit r/ — )• vr" or A — )• —1"*" will be 
discussed at length in Section [5l (b) The entanglement entropy of two spins m+i} at the 
Free Fermion point, where the correlation functions (see e.g. j43j for a derivation) are given by: 

rmi r m+l -i 

l^z. \ _ 2((-i)'"-i) , , + . _ i-ir n r(fc)^ A r(fc)^ 
(fTi^T^+i) - ^2^2 and {a^ a^^^) - [[^ T{k-\)T{k+\) r(fc-i)r(fc+i) ' 



Finally, we note that the correction terms for the asymptotic behaviour of the entanglement 
entropy at large m and generic values of the anisotropy parameter can be obtained by employing 
the results of [S]. In this paper, the asymptotic behaviour of the two-point functions (cfo",^^]^) 
and (cfcr^^^) was obtained using perturbed CFT techniques (the results of [2] were even more 
general, as they considered time-dependent correlation functions). 
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Using 



r 



A 



V 



2(i-f)M2vTr^ 



exp 



dt 







sinh(-t) 



t ysinht cosh(l — tt 



(3.59) 



we obtain the following leading asymptotic expression for the entropy: 

6t) ') ^""11 



•^ll.m+l} ~ 21og(2) 



m T 8m ^ 



4^ +0{m ^' ,m ^ 'f,???, '^), (3.60) 



for T] G (0, tt) 



4 Quantum states for chains of infinite length 

Infinite-length chains can be described using Hilbert spaces. However, contrarily to the case of 
finite-length chains, the infinite-dimensional Hilbert spaces Ha occurring in the limit N ^ oo 
depend, in a sense, on the parameter A characterising the Hamiltonian ()3.36|! ll. More precisely. 
Ha is formed by the completion of the set of finite-energy eigenvectors of the Hamiltonian Ha , 
and these are sets of very different vectors for different values of A. In particular, the Hilbert 
spaces T-La, T~La' are orthogonal for A ^ A'; e.g. in the limit of infinite volume, the overlap 
a('0|'0)a'/ ■\/a(V'I^)aa' (V'I'0)a' between the two ground states tend to zero. Due to this, the 
description of infinite-length quantum chains via Hilbert spaces is not very well adapted to the 
study of quantum averages of localised observables. Indeed, the quantum average of an operator 
that factorises to the identity on all sites except finitely many is well defined on the Hilbert 
spaces T-La for any A (in the range — 1 < A < 1 considered), and smooth as a function of 
A. The Hilbert space description does not provide any natural topology in agreement with 
continuity in A. The study of the limit A — t- — 1"*" requires a more adapted description. 

For Hilbert spaces of finite dimension, knowing all quantum averages uniquely fix the quan- 
tum state. From this, one way of describing quantum states (pure or mixed) that is more 
convenient when considering the infinite- volume limit, is through linear functionals representing 
quantum averages. Let us consider the infinite set of vector spaces {Vi : i G Z}: this is the 
set of sites of the chain as usual, except that the index can be both positive and negative (i.e. 
there is no boundary). For any finite block of sites A, let us consider the complex vector space 
£a of operators supported on A: the space (8'jg^End(Vi) with action like the identity on all 
sites j A. Let us further consider the complex vector space S = Ua^a of finitely-supported 
operators; explicitly, these are the finite linear combinations of finite products of elementary 
matrices . 

Definition 4.1 The space of quantum states on the infinite-length chain, which we will denote 
by J-, is the space of all linear functionals tp : £ ^ C that are real-linear: 

^{Ef = ^{E^) {E G £) (4.61) 

and normalised 

ip{l) = 1. (4.62) 
''Of course, as a Hilbert space, the space with an infinite countable basis is unique. 
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The physical meanmg of ip{E), for hermitian is simply the average of the observable E in 
the physical state represented by "ijj. 

For infinite length, this definition of a quantum state is extremely useful, because it puts in a 
common space the quantum states of any quantum chain with well-defined correlation functions: 
for instance, the XXZ chains for —1 < A < 1. Indeed, to every ray {ajV') : a G C} lying in the 
Hilbert space H.^ for some A, we associate a quantum state in T defined by 

This definition of the space of quantum states is also very natural from the viewpoint of the 
Renyi entropies S^°°^'(n) for integer n. Given a quantum state ijj ^ J-, we can construct the 
linear functional ^ := tp^"^ on S^"^, and the Renyi entropy of ip with respect to a finite set of 
sites A is expressed using our replica permutation operators: 

Sf-y\n) = log (Ta)) . (4.64) 
1 — n 

Note that indeed, the permutation operator Ta associated with a finite block of sites ^ is a 

finitely-supported operator. Naturally, since the density matrix is completely determined by 
the fundamental blocks on the sites A, it can also be defined in the infinite-volume case. In 
particular, we have 

Tv{p-^) = ^{rA). (4.65) 

Recall the concept of local-unitary transformations introduced in section [2l This concept 
can naturally be included into the infinite-length setting. Given disjoint blocks of sites Ak, 
k = — oo, . . . , oo, a block-defined transformation is a doubly-infinite sequence U = (Uk : k £ Z) 
of elements Uk € "fA^ that are invertible, V k. To any such U we can associate an adjoint 

action on £, given by Ad U{E) := Hfc U^^EUk, E € £. The product is finite, hence this is a well- 
defined action. This gives rise to an action on quantum states, defined hy U ■ ip = ip o Ad U for 
any tp £ T. K local-unitary transformation is a block-defined transformation, where Ak = {fc} 
all Uk are unitary. It stays true in the infinite-length setting that Ta is local-unitarily invariant. 

In order to discuss convergence of quantum states in J^, we define a topology (in fact, 
a geometry), via compact convergence of fundamental blocks. This topology in fact works 
(and will be used) more generally for linear functionals without the normalisation condition 
(^(1) = 1). This is inspired by the usual compact convergence topology in function theory. 
Consider a sequence of functions on T parametrised by an integer A; > 0, defined by 



4(V'i,'(A2 



\ {ei,€[:t=-k,...,k} 



\i=—k / \i=—k / 



2 



This is a natural definition, because it is invariant under a change of orthonormal basis in the 
space of linear operators on the subchain of 2 A; -|- 1 sites centered at (with the inner product 
given by {E,E') = Ti{E'^ E')). Consider then the distance function 

I \ 2-^4('0i,'02) 

Definition 4.2 The compact convergence topology on the space of quantum states T is that of 
open balls induced from the distance function (f^.gg[ ). 
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This topology is the most convenient when considering finitely-supported observables. In- 
deed, by the usual arguments of compact convergence, one can see that every Cauchy sequence 
{ipj G J-" : j = 1, 2, . . .) gives rise, for any finitely-supported operator E, to a Cauchy sequence 
(ipj^E) : j = 1, 2, . . .). In other words, the limit of a Cauchy sequence in is in J^: the space of 
quantum states is complete. Moreover, again by standard arguments of compact convergence, 
if, for a sequence (■0^- £ T : j = 1,2, . . .), all fundamental blocks tpj (jXi^mi -^i^'*^*) converge as 
J — 7- oo, then the sequence converges, and the limit is a quantum state whose action on ele- 
mentary matrices agrees with the limit of the fundamental blocks. That is, this is the topology 
that associate convergence of quantum states to convergence of averages of all finitely-extended 
observables. 

We remark that this topology also allows us to describe the infinite-length limit of quantum 
chains. Indeed, to any sequence of vectors lip)^^^ in the space spanned by the sites from —k to k, 
for /c = 1, 2, 3, 4, . . . , we can associate a sequence V'fc of elements in by V'/fcC-^) = ^'""^(V'l-^'IV')^''"^ 
if E € £{-k,...,k}y and ipki^) = otherwise. If \ip)^''^ describe ground states, and if the infinite- 
volume limit of the chain exists, then this sequence converges, hence define a state limfc_^oo ^fc ^ 
T. The same procedure can be applied for any sequence of linear functionals on spaces Sa^. with 
growing such that lim^^oo = ^- Let us also note that in this context, there is no a priori 
clear distinction between pure states and mixed states in infinite volume. Indeed, for any density 
matrix p on a finite number of sites, it is always possible to construct a pure state on twice as 
many sites, such that the reduced density matrix equals p. Hence, to a sequence of density 
matrices leading to a converging sequence of quantum states, we can associate a sequence of 
vectors giving rise to the same limit quantum state (and vice versa). 



5 The limit A ^ -1+ 

We now wish to study the state occurring in the limit A — )■ — 1^ of the infinite-length XXZ 
model in compact convergence topology. Using the calculation Rules [2Tl as well as exact results 
for fundamental blocks of the XXZ models found in the literature, we provide a full description 
of the reduced density matrix associated to any finite subset of sites. Our main result is that 
the reduced density matrix associated with a number m of sites (no matter what their positions 
are) has only m + 1 non- vanishing eigenvalues which are given by 

1 / TTi \ 

Afc = — ( ^ 1 with k = 0,...,m. (5.67) 

Therefore, for m odd each eigenvalue is twice degenerate, whereas for m even all eigenvalues but 
one are twice degenerate. 



5.1 The XXZ quantum state in the Hmit A — — l"*" 

From the explicit formulas obtained in \39\ \3E[ [371 140j . it follows that the limit A — )■ —1^ of 
the fundamental blocks associated with the ground state of the Hamiltonian ()3.36p exists and is 
finite. More precisely, specializing the multiple integral formulae obtained in |39 [ 1381 [371 140j . we 
find that in the limit A — t- — l"*", the non- vanishing fundamental blocks associated to an index 
set A = {ji,j2, ■ ■ ■ ,jm} (with jk / jk' for k / k') reduce to 

lim {E'''''Ef''^ ■ ■ ■ E'"'''^) = — TT (-1)^ (5.68) 
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where B is the subset of sites at which either an operator E or an operator E sits: 

B = {h : efc / 4} C A. 

Note that the indices ji, ■ ■ ■ ,jm in formula (j5.68p are not necessarily consecutive, contrarily to 
the formula for fundamental blocks in \37\ [36] . It is simple to check that ()5.68p is consistent 
with the relation 1 = Ej^ + Ef. 

Let us consider the ground state |V')a of the Hamiltonian ff^- For ^ £ ("Ijlli the corre- 
sponding quantum states all lie in T. Moreover, by formula (j5.68p . all fundamental blocks 
converge in the limit where A tends to —1 from above. Hence, the limit lim/'^_j._i+ -i/'^ exists 
and is a quantum state of the infinite-volume chain. We will denote this limit by 

V'. := lim G F. (5.69) 

This quantum state is completely characterised by formula (|5.68p : 

( E'-'^'E'^''^ ■ ■ ■ E""''A = 

We would now like to study this state by computing its Renyi and von Neumann entangle- 
ment entropies. As is clear from (j5.70p . a very special feature of this state is that its correlation 
functions are extremely simple. Hence, we will be able to find analytic expressions for the Renyi 
entropies and deduce the exact eigenvalues and degeneracies of the reduced density matrices; 
these are inaccessible for generic values of A. 



— TT(— 1)"' equal numbers of E^"^ and E'^^ 
otherwise. 



5.2 Renyi and von Neumann entanglement entropies of the state 

The Renyi entropy ()1.3p can be evaluated using the replica permutation operators via (|4.64p . 
The present subsection provides a proof of the following proposition, which is in agreement with 
the first point stated at the beginning of this section, equation ()5.67p . 

Proposition 1 Consider the state obtained in the limit A — )■ —1^ of the ground state 
in the compact convergence topology. Its Renyi and von Neumann entanglement entropies 
associated with any set of m sites are given by 

and 

5.^™lo,2-iE(»)l„.(7). (5.72) 

fc=0 ^ ^ ^ ^ 

The large-m asymptotics of these quantities are given by 

SrV) 4 |iM_ . O 5. 4 lo, ( = ). i . O (,„-) . (.73) 

The leading term in the asymptotics of the Renyi entropy is therefore independent of n. 
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Let us start with the simplest non-trivial case n = 2. The specialisation of (j2.18p to 7t, = 2 
immediately gives the formula 



e',ei = l 

Notice that many of the correlation functions in this sum are vanishing: a fundamental block is 
non-zero if and only if it involves equal numbers of matrices E^^ and E'^^. 

Given the structure of the product (|5.74p it is easy to see that the products of (—I)-' in ()5.68p 
will not play a role, as matrices E^"^ and E'^^ appear at exactly the same sites in both correlation 
functions, so that all phase factors cancel out. Therefore, the sum (|5.74p reduces to 

^ / 9 X number of non-vanishing m-site fundamental blocks , „ , 

(pa) = ^ • (5.75) 

The number of non-vanishing m-site blocks can be easily obtained by simple combinatoric argu- 
ments. A generic non- vanishing fundamental block will have s operators E^^, q pairs of operators 
E^"^ and E'^^, and m — s — 2q operators E'^'^ at the remaining sites. The total number of such 
blocks is simply 

(5.76) 



therefore 



s\iq\y{m- s-2q)\' 

[t] m-2q . -n n I ™\ 

22™ ^ ^ s!(g!)2(m-s-2g)! ^T{l + m)' ^ ' 

The Renyi entropy ()1.3p is given by 

(5.78) 



gr^X2) = -log( rrn^K 

\ V'^r(l + m) 



and, for m large behaves asymptotically as 

^m ""H2) = I log(vrm) + ©(m^i). (5.79) 

The calculation can be extended to higher values of n. The Renyi entropy for generic n 
will be obtained form a sum over an n-fold product of m-site fundamental blocks. A feature of 
Rules [2T] is that for every site, there must be as many matrices E^'^ as there are matrices E'^^ 
distributed amongst the n copies. Hence, as in the n = 2 case, the (—1)-' factors in (|5.68p all 
cancel out when multiplying the n fundamental blocks together, because any possible factor — 1 
will always be present an even number of times. This means that once more, the computation 
reduces to combinatorially working out the number of non-vanishing products of fundamental 
blocks that occur. 

A simplification to this combinatorial calculation is obtained from two general consequences 
of Rules 12.11 First, once a choice of elementary matrices has been made for the first n — 1 
fundamental blocks, the choice for the n**^ one exists and is unique. This is because the rules 
say if, on any given site, we have the elementary matrices E''+^+ at copy a + 1 and at 
copy a — 1, then the elementary matrix at copy a is fixed to be . Second, if the first n — 1 

fundamental blocks are non- vanishing and satisfy the rules, then the unique choice of the last 
block also is non- vanishing. This is because there is the same number of matrices E^'^ as there 
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is of E'^^ in total in the first n — 1 blocks (since this is true in every block). Since at every site 
there can only be, in the n — 1 first copies, an equal number of E^'^ and E'^^ or a surplus by one 
of either E^"^ or E'^^, then the number of sites where there is a surplus of E^"^ must be the same 
as that where there is a surplus of E"^^. Hence, Rules ETT] imply that on the last block, there 
will be the same numbers of E^'^ and E'^^. Therefore, putting these three observations together, 
in order to count the number of non-vanishing products of n fundamental blocks occurring, we 
only need to count the number of choices of elementary matrices for the first n — 1 fundamental 
blocks that make these blocks non- vanishing and that satisfy Rules 12.11 

Let us now analyse the case n = 3, where the quantity Tr (p^) is given by a sum over triple 
products of m-site fundamental blocks. We can depict Rules [TT] through diagrams such as 

1st fundamental block Ej^ Ej^ 

v/ \ v/ \ 

2nd fundamental block Ej^ or Ej^ Ej^ or Ej^ 

4^ ^ ^ ^ 

3rd fundamental block Ej^ or Ej'^ Ej^ or Ej'^ 

and similarly for Ej"^ and E^^ by exchanging super-indices 1,2 everywhere. 

Let us consider the first fundamental block. Let it contain s matrices E^^, q pairs of matrices 
^12^ ^21^ g^j^^ g_2q matrices E"^"^ . Each such fundamental block can appear in conjunction 
with various choices of the second fundamental block. According to the rules above, every 
matrix £'^^ or E^"^ in the first can be followed by a matrix E^^ or E'^^ in the second. Similarly, 
every matrix E"^"^ or i?^^ in the first can be followed by a matrix E"^"^ or E^"^ in the second. The 
number of second fundamental blocks containing a pairs of matrices i?^^ and E"^^ is given by 



s + q \ I m — s 
a I \ a 



(5.80) 



Summing over all possible values of a gives 

min(s+5,m— s— (jf) 



E 



s + q \ I m — s — q \ _ I m 
a j \ a I \ s + ( 



(5.81) 



a=0 

Therefore, the total number of distinct triple products leads to 

[^] m-2q , 

\Pa) - ^ Z> s!(g!)2(m-s-2g)! \s^q 



2'"0iT(l -Km) 

The leading asymptotic behaviour of the Renyi entropy can be extracted by exploiting the fact 
that for large m 



3F2 - m,-m, -y| , |l, ^ - m| ; 1^ ~ 3F2 (^|-m, -m, -y | , {1, -m} ; 1^ 

«m m^ , ^ \ ml , , 
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Substituting this into (j5.82|) and employing Stirling's approximation for the Gamma functions 
we obtain, 

5^'"^'(3) = \ log(^) + 0(1, m-i). (5.84) 

Unfortunately, the approximation above is too crude and does not allow to extract the correct 
constant 1/4 log (3) as predicted by ()5.73p . This is the same leading behaviour as in (j5.79p which 
is incompatible with the CFT prediction according to which the coefficient of log(m) should be 
a function of n. 

It is possible to extend the previous combinatorial arguments to higher (generic) values of 
n. Again, from Rules I2.H we know that it is from the set of matrices E^^ and E"^^ in the 
fundamental block of copy a that emerge all matrices E^^ and E^"^ in the copy a + 1. The same 
is true if all super-indices 1,2 are exchanged. Let us consider the a^^ fundamental block, for 
a = 1, 2, . . . , n — 1. Let it contain Sa matrices E^^, pairs of matrices E^'^,E'^^, and m— Sa — '^Qa 
matrices E'^'^ (with s := si and q := si). Since Qa+i matrices E'^^ are chosen from the Sa + Qa 
matrices E^^ and E^'^, the rest must be all matrices E^^, so we must have 

Sa+l = Sa + Qa - Qa+l- (5.85) 

Hence, we find 

Sa+l + Qa+l = Sa + Qa =^ Sa + qa = S + q^ O. (5.86) 

Therefore, from expression ()5.80p . the number of fundamental blocks at copy a is given by 

s + q\f m-s-q\ 
qa J 

Summing over all possible values of qa for a = 2, . . . , n — 1 and using formula (j5.8ip gives 

^ [t1 m-2q ^1 , s n-2 

^(p1) = ^EE sliqiy{m-s-2q)\[s7q ) ' ^^'^^^ 

Prom the definition (jl.Sp it is clear that the formula above should give exactly the sum over 
all n-powers of the eigenvalues of the density matrix associated to the state . However, this 
is not obvious from ()5.88p in its present form. In order to make this structure clearer it is useful 
to rewrite (|5.88p in a slightly different way. We will introduce the new variable k = s + q and 
write the binomial coefficient to the power n — 2 as a product of powers n and —2. The sum 
becomes, 

f— 1 

yt'A) 2^"^ ^ ^ (k-q)l(ql)^(m-k-q)\ml\ k J " ^ ^ ^ 

q=0 k=q J./ \i / \ J/ \ / 

Exchanging the two sums we find the equivalent expression: 

^ m m—k TM\2 / \" i m. / 

T r " ^ = — V V [k\(m-k)\Y ( ^ \ = J_ V ( 1 r5 qn^ 

^^{Pa) 2-^'^^^{k-q)\{qiy{m-k-q)\ml\k J 2^'^^\k ) ' 

where the last equality follows from the remarkable identity 

m—k 



{k\{m-k)\f 



(k — q)l{ql)'^{m — k — q)\m\ 



(5.91) 
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The result (j5.90p clearly identifies the non- vanishing eigenvalues of the density matrix as (j5.67|) 
and provides a much simpler expression for the von Neumann and Renyi entropies, given by 
(I5JT]1 and ([5772]) . 

Even though we have not been able to carry out the sums explicitly, it is possible to extract 
the leading behaviour of (j5.7ip and (j5.72p for large m. First, we transform the last sum in (|5.90p 
into a multiple integral using a generating function procedure: 

1 1 r""^ . / 



n— 1 n—1 



^ ^ ^ i=i ^ 



cos 



(5.92) 



Then, we evaluate the large-m asymptotic using a saddle-point analysis. The main contribution 
to the integral at large-m will be obtained when the integrand is evaluated at ~ for all j. 
The leading asymptotics is evaluated by using cos(^/2) = ^-^"^ I'^^Oie^).^ 



1 



(27r)"-i 
1 
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JJ^ dQj exp 



(2vr)' 



1 



n-l 



n-1 



d^j exp 



n-l 

i=i Vi=i 



(5.93) 



i=2 



n-l 

2 \~ 



(5.94) 

/n \7rm / 

Corrections coming from the 0(^'^) terms in the exponential give rise to a factor 1 + 0(m~^). 



5.3 Discussion 

In order to discuss the results of Proposition [H we construct the state ^_ in a different way. 
Let be the operation, on the space of finitely-supported operators, that projects onto those 
keeping the z component of the spin: 



£ ^ £ (linearly) 



J' J-x 



Further, let U be the following local-unitary operator: 



U 



n ' 

je2Z+i 



and let be the following quantum state: 



where V'ea:;^ • ^{j} — ^ C is the average, on site j, obtained from the vector {ip^ 



(5.95) 



(5.96) 



(5.97) 



t>i + I i 

^ representing a spin pointing in the positive x direction (this is the eigenvector of 
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with eigenvalue 1). Note that the right-hand side of (j5.97p is a converging doubly-infinite tensor 
product in the compact convergence topology. Then, it is a simple matter to see that 



(5.98) 



Indeed, from (jOTl) we notice that ^p^^.j i^Ej'^ j = 1/2 for all e,e'. Using that Ad U \^Ej'^ 

{—lyEj''' for e ^ e', and that Ad U (^Ej''' ^ = Ej'^ for e = e', we immediately find that the right- 
hand side of ()5.98p reproduces formula ()5.70p . Note that in ()5.98p . the order of the operations 
U ■ and ¥z is irrelevant, since they commute (i.e. [Ad U , P^] = 0). 

The form ()5.98p of the state points to the following construction. For each A: = 0, 1, 2, . . ., 
consider the Hilbert space 0^^_j^Vi, and the vectors where all spins point in the positive x direc- 
tion: iV'e^;)''^'^'''"'^^ •= HiL-fc IV'ej;)*- Of course, these vectors can be written as linear combinations 
of vectors with spins pointing in z directions (here, we use the standard values ±1/2 =t, i for 
the spin variables Sj, instead of the numbers 1,2 respectively): 



v(2fc+l) 



1 



iV2) 



2k+l 



{Si = ±l 



.,fc}' 



Now consider the Hilbert space (^i=_kVi(^Vcnv with Vcnv = C^''+^. The space Vc 



of as a spin- ( A; -|- 1/2) module, with orthonormal basis vectors \S)q 
of spins 5/2. Take the vectors \^pj^) as follows: 



2k+l 



E 



Si) 



-k,...,k}'' 



n Sg{- 



E 



env can be though 
A:-l/2,...,A; + l/2} 



(5.99) 



Essentially, we have adjoined to every vector in the sum an additional degree of freedom with 
a spin exactly opposite to the total spin of the original chain (the sites from —k to k), in such 
a way that all vectors have total spin (in the z direction). In the sub-sum of vectors where 
the original chain has a given total spin, the adjoined vector just factorises out, so that when 
calculating averages, the resulting sub-sum of matrix elements is unchanged. Hence, the average 

in [ipi^) of any product nj=-fc ^Z' ^ tha-t preserve the total spin of the original chain is exactly 
equal to the average in I'fpex)^'^^^^^ ■ However, the average of such elementary matrices in {tpj^) is 
zero if they do not preserve the total spin. 

As usual, from \ipf.) we can form a linear functional ipf^ on £{-k,...,k}- By the discussion above, 
we have 



n4 



(^''+^)(V'eJ JJ Ej"^'\'iljgJ^'^''^^'> total z component preserved 

j=-k 

otherwise. 



We may now take the limit 



lim tp^. 



(5.100) 



Since the averages in the states lip e^)^'^^^^^ stable as k increases, and since the linear func- 
tionals act on increasing subsets of the chain whose limit is the whole chain, we find that the 



limit exists in the compact convergence topology. We have 



ijj^^ P^, hence 



(5.101) 
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The state ipg^^ is factorisable, hence has entanglement entropy equal to 0. The higher entan- 
glement entropy of ()5.98p (in particular, the diverging behaviour at large m) comes from the 
interplay between this factorised state where spins point in x directions, and the condition of 
having a total z-component equal to 0. 

Construction (|5.99p . (|5.100p . (jS.lOip gives some insight into this interplay. It makes it clear 
that the extra entanglement entropy generated by the presence of the operator P^, i.e. by the 
condition of preserving the z-component, is due to an extra entanglement with the additional 
degree of freedom Vcnv This degree of freedom can be interpreted as an "environment", which 
couples with the "system", the sites from —k to k, via the z-component of the total spin. The 
system can be seen as a random collection of spins up and down, all configurations having 
equal probabilities. Such a state has zero entanglement entropy with respect to any subsystem. 
However, the coupling generates entanglement between the system and the environment, which 
is picked up when measuring the entanglement between a subsystem and the rest. Since the 
environment is infinite-dimensional in the infinite-A; limit, its capacity for entanglement is infinite, 
hence as the number of sites m of the subsystem increases, the entanglement entropy increases. 
On the other hand, since the environment does not couple to each individual spin separately, 
the entanglement entropy increases at a rate that is less than linear in m. The number of blocks 
of system states that couple to separate environment states is the number of possible values 
of the z-component of the total spin. For m sites, this is m + 1, hence we can indeed expect 
a logarithmic behaviour in m for the entanglement entropy at large m. Since the blocks have 
different sizes, they have different probabilities of occurring, hence we may indeed expect b log m 
for some b < 1. In fact, the size of a block with total z-spin equal to k — m/2, for k between 

^ 1 . This explains the eigenvalues of the reduced 

density matrix, from which we derived the (1/2) log m behaviour. 

Besides explaining the features of the entanglement entropy of the state ip_, this construction 
also gives an insight into the way the ground state of the XXZ model tends to as 
A —7- — l"*". First, note that U ■ ip^^ is a ground state of the XXZ model at A = —1, because 
(as is well known) the Hamiltonian of this model is simply obtained from a U transform of 
the Hamiltonian of the ferromagnetic XXX model - the ground state energy of the latter is 
highly degenerate, and ipg_^ is a possible ground state. We expect that as A — )• — 1+, there is 
an increasing characteristic length ^ such that: 1) on scales much below ^, sites are essentially 
randomly up and down (in the z-direction) , with appropriate coefficients so that they locally 
reproduce the ground state U ■ ip^^ of the XXZ model at A = —1; and 2) on scales above ^, the 
total spin is more likely to be zero, as it should be for the XXZ model for — 1 < A < 1. As 
A — )• —1"'", the "inner" sites, on any finite scales, are described by the state U ■ ipe^, and the 
"outer" sites, beyond ^ — )• oo, are collectively described by the single degree of freedom V^nv, 
whose entanglement with the inner sites, necessary to keep the total spin to zero, projects onto 
spin-preserving operators. 

An immediate feature of large-m scaling of the Renyi entropy of the state "0- is that it is 
in disagreement with the predictions coming from conformal invariance of critical points. For 
instance, for a region A oi m consecutive sites, CFT predicts 

5^-yi(n) ^(1 + n-^) log m + 0(1) 

where c is the central charge. This is different from ()5.73p . which is independent of n. Yet, the 
scaling (|5.73p is logarithmic, hence there is a similar type of scale invariance (or covariance): 
a region twice bigger sees its entanglement entropy added by half a "unit", 1/2 log (2). The 
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scale ^ should give rise to a new universal function interpolating between this non-conformal 
scale-invariant behaviour, and the conformal scaling. That is, it should be possible to define a 
new scaling limit, looking at the entanglement entropy for a region of length m that scales like 
m = while the limit A — >■ — 1+ is taken. The result is likely to be /(a) logm, for a function 
/(a) interpolating between /(O) = 1/2 and the coefficient predicted by CFT, e.g. /(oo) = c/3 
for the von Neumann entanglement entropy of a connected region, with the central charge c = 1. 

Another important feature of the entanglement entropy of the state V'- is that it does not 
depend on the actual positions of the sites of the region A. As will be clear from our future 
work [29], this is a feature of any ground state of the XXZ model at A = — 1. It turns out that 
the state V'- is just one of these ground states. The field theory interpretation is that the twist 
field has dimension (corresponding to a central charge 0), but that there is an entanglement 
entropy contribution coming from the degeneracy. 

6 Conclusions and outlook 

In this paper we have described in detail a new approach to the computation of the bi-partite 
entanglement entropy of quantum spin chain models. Our approach employs the "replica trick" 
widely used in the CFT context \\b\ I16j and also, more recently, for integrable QFTs [21]. This 
approach consists of considering a "replica" version of the theory for which one wants to compute 
the entropy, that is, a new theory consisting of n non-interacting copies of the original model. 
It turns out that computations of the entropy are often more feasible in the replica model than 
in the original seemingly simpler theory. The reason for this is the existence of new symmetries 
of the replica model which enable extra fields (in the QFT case) or extra local operators (in the 
quantum spin chain case) to exist in the Ji-copy model. These fields (operators) are related to 
the symmetry of the replica theory under cyclic permutations of the n-copies. 

For integrable QFTs the fields referred to above were introduced in [24j and named branch 
point twist fields. It was shown that the entanglement entropy can be expressed in terms of 
their two-point functions which lead to various computation of the entropy which have been 
summarized in the review |28j . 

In this paper we have identified the basic operators in terms of which the spin chain equivalent 
of the QFT twist field can be expressed: they are local cyclic replica permutation operators. Each 
operator % acts on one site of the chain i by cyclicly permuting the spins of the n-copies at that 
particular site. Ground state correlation functions of products of such operators over different 
sites of the chain play the same role as two (or higher)-point functions of branch point twist 
fields in QFT. 

We have described a representation of the replica permutation operators for quantum spin 
chains of general spin in terms of fundamental 2x2 matrices E"^- acting at site i, copy a of 
the chain and proved that such representation does reproduce the cyclic permutation action on 
quantum states described above. We have also shown that the permutation operators satisfy 
exchange relations with respect to other local operators of the chain, in much the same way as 
the twist field in QFT. In fact, as in QFT, these exchange relations can be used as starting point 
for the construction of the permutation operator in quantum spin chains. 

We have demonstrated the working of our approach by evaluating the bi-partite entanglement 
entropy of one and two spins in an infinitely long spin-^ quantum spin chain characterized by 
an antiferromagnetic ground state. We find that the computation of the entropy is recast into 
combinatorial problems, which become less tractable as the number of spins is increased. It 
will be very interesting to address this problem in the future and to establish whether or not 
our approach may be useful for the study of the entanglement entropy of large subsystems. 
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From our approach it is clear that the larger the subsystem the more operators will be involved 
in the correlation functions which enter the expressions of the entropy, so that advancement 
in the understanding of the large-distance behaviour of such correlation functions (see [55] for 
recent progress) goes hand in hand with the success of our approach to evaluate the large size 
asymptotics of the entropy. 

Despite the difficulties emphasized above we have been able to identify a particular quantum 
state for which both the von Neumann and Renyi entropies are non-trivial but can still be 
explicitly computed for any subsystem size and any values of n. This quantum state, which we 
denoted by is defined by its correlation functions which are obtained as the A — )• — 1"*" limit 
of those of the spin-^ XXZ quantum spin chain in the antiferromagnetic regime |39 [ \'6H \ [57 1 14U]. 
Our analysis has allowed us to identify all the eigenvalues of the density matrix associated to 
this quantum state and to find the exact large subsystem asymptotics of both the von Neumann 
and Renyi entropies. In both cases the leading term scales with the logarithm of the size of the 
subsystem as for critical systems. However the coefficient of this logarithmic term is independent 
of n for the Renyi entropy in direct contrast with the well-known dependence found for CFTs. 

It is tempting to speculate that this logarithmic but non-conformal behaviour may be some- 
how related to the unusual behaviour at A = — 1 encountered in [5]. In this recent work it is 
argued that such unusual features are characteristic of essential singularities. We must however 
stress that the quantity evaluated in [S] and in the present work are rather different in nature, 
as well as the characterization of the quantum state of the chain, so that the unusual features 
observed in both cases are not a priori related in an obvious way. Yet, it is possible that the 
description that we gave of how the ground state approaches the A = — 1 point from the region 
A > — 1 may be of use in understanding the results of [5]. For instance, it may be that the 
characteristic length ^ could be replaced by the correlation length in describing the approach of 
the point A = —1 in [5]. The analysis of the limit A — )• — 1, and the results of [5], naturally beg 
the question as to the entanglement entropies of the infinitely-many ground states of the XXZ 
model at A = —1. As we said, it turns out that the state ^p_ is just one of these ground states, 
although a very particular one, highly entangled. Likewise, it is possible to construct a ground 
state for every possible finite asymptotic value of the entanglement entropy and it may be that 
behaviours observed in [5] result from approaching these various ground states. We will come 
back to the study of the entanglement entropies of A = — 1 ground states in a later work |29j . 
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